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Abstract. Bounds on the trace mappings defined on the Sobolev space W^(Q) 
and the space LD(Q) of integrable stains are obtained. Such bounds corre- 
spond to stress concentration — the ratio between the maximal stress in a body 
and the maximum of the traction applied to its boundary. The analysis lead- 
ing to the bounds may be described in the mechanical context of stress theory 
and stress concentration. 



1. Introduction 

This work considers some mathematical aspects of stress concentration. Stress 
concentration is a term used by engineers to indicate the increase of stresses from 
some expected values due to deviation of the geometry of the body from an ideal- 
ized simple one. In practice, stress concentration factors are used by engineers to 
indicate the ratio between the maximum of the stress for the actual body under 
consideration and the maximal stresses calculated for the simplified geometry for 
which simple formulae of strength-of-materials are used traditionally. Stress con- 
centration factors are usually compiled for homogeneous, linear, isotropic, elastic 
solids for various typical geometries (see for example [7j)- Their values are obtained 
by solving the equations of elasticity analytically, by numerical approximations and 
by experimental methods. 

In recent work, [HI Eli we formalized the notion of stress concentration mathemat- 
ically and generalized it. The idea in formalizing the stress concentration factor is 
to regard the stresses for the simplified geometry as traction boundary conditions. 
For the example of a bar under tension F, the "nominal" stresses F/A, where A 
is the cross section area, are regarded as boundary conditions at the ends of the 
bar in agreement with the engineering notion. Thus, for a given surface traction 
distribution t on the boundary of the body fl and a stress tensor a in equilibrium 
with t, the stress concentration factor is defined by 

ess sup.,. \a(x)\ 
ess sup,, \t(y)\ 

Here, |<r(a;)| is a the norm of the value of the stress at x. To evaluate it we use 
some norm on the vector space of matrices. The value of Kt,a clearly depends on 
the norm used and various norms are discussed in Section 4. Using the essential 
supremum, we ignore high stresses on sets of zero volume. 
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Noting that without specifying particular constitutive relations there is a class 
St of stress fields that are in equilibrium with t, we define the optimal stress 
concentration factor as 

K t = inf K t , a . (1.2) 

Next, realizing that an engineer usually does not know a-priori the nature of the 
loads acting on a body exactly, the generalized stress concentration factor is defined 
as 

v v f t ess sup x£f2 \a{x) \ \ 

A=supA t = sup< mi . . > (1.3) 

t t [vest ess sup yega \t(y)\ J 

where the supremum is taken over all traction fields — essentially bounded vector 
fields on dQ — i.e., over all t £ L ca (dQ 1 R 3 ). It is noted that A' is a purely geometric 
property of the body J7. 

In we related the generalized stress concentration factor to the norms 
of the trace mappings on Sobolev spaces and LZ?-spaccs. It turns out that for 
the formulation of equilibrium and stress theory, particularly in the context of 
stress concentration, the Sobolev space Wi(f2,M. 3 ) and the related LD(Q) space 
are especially useful. We recall that the space LD(Q) contains integrable vector 
fields w such that the components of their associated stretching, or (infinitesimal, 
linear) strain, 

e{w) = i(Vw + (V W ) T ), e(w)ij = ^(w itj + wjj, (1.4) 

are also integrable (see [T21ITUllTn i31). 

Assuming that J? is an open subset of K n having a C 2 -boundary, for both spaces 
one has a well defined, linear, bounded trace mapping 7, such that for every vector 
field w defined on i?, j(w) is a vector field defined on df2 satisfying the following 
compatibility condition. For every continuous vector field u defined on the closure 
f2, 7 acts as the restriction to the boundary, i.e., 

j(u\a) = u\an- (1.5) 

In [5] we have shown that if we ignore the requirement that the total force and 
total torque on every subbody of fl vanish, then, 

#=11711, for j:Wl(n,M. 3 ) — ► l}(dn, M 3 ). (1.6) 

For the case where the total forces and torques on the various subbodies do vanish, 
a more detailed analysis is required (see [Hj)- Letting be the finite dimensional 
vector space of rigid vector fields, one has to consider the quotient spaces LD{fl)/.S# 
and L 1 (9i7, R 3 )/^. For these spaces, one can define an induced trace mapping 

7/^: LD{Q)/M — ► L l {dn,W)/M, (1.7) 

and it turns out that 

K = h/M\\. (1.8) 
Thus, apart from the mathematical interest in estimates on ||7|| and Ht/^H, such 
estimates are very significant in stress analysis. It is our objective here to estimate 
these constants. Our method of estimation has a mechanical flavor. In a way, it is 
dual to the analysis leading to the relation between generalized stress concentration 
factors and the norms of the trace mappings. 

For the particular case of the Sobolev space Wl(f2), M. Motron obtained recently 
jSj some estimates on the bounds. The method used here is different and is based 
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on the maximum principle for the Dirichlet problem. The bounds we obtain give 
a concrete estimate for ||7||. Subsequently, we extend our method to the space of 
LD-fields. Specifically, recalling that the LZ?-norm is given by 

\\w\\lD = ^2\\wi\\ L i + Y)\ S ( W )im\\L^, (1-9) 

i i,m 

we obtain bounds on the constants A and B such that 

J\ 7 (w)\^Aj\w\+Bj\e(w)\ (1.10) 

an n n 

for all LD-fields w (so that ma,x{A, B} bounds ||7||). 

Section 2 presents the results and methods of [S] relating stress concentration 
to the norm of the trace mapping. The last subsection discusses the simplification 
to the case of sourceless vector fields (rather then stress tensors) which may serve 
as motivation for studying the norm of the trace mapping for the Sobolev space 
Wl(O). In Section 3 we introduce the basic method for obtaining the bounds on 
the trace mapping for the Sobolev space using harmonic vector fields. The main 
result of this section is Theorem 13.41 As background material for the discussion 
of the bounds for the trace mapping on LD(f2), Section 4 presents some standard 
results on norms of matrices. These are significant in the mechanical context as 
they are used on the space of stress matrices. For example, a yield criterion is 
usually a seminorm on the space of stress matrices. Section 5 presents additional 
preparatory material — the optimal boundary values for the stresses for a given 
boundary traction field. Section 6 studies the bounds on the norm of the trace 
mapping for LD(fi) using the method of harmonic tensor fields and the central 
result is Theorem 16.11 Finally, the concluding remarks of Section 7 discuss the 
mechanical interpretation of the preceding analysis. 



2. Generalized Stress Concentration and 
the Norm of the Trace Mapping 

2.1. Basic definitions and notation. We consider an open set Q C M 3 , where 
in some sections, the presentation is in the setting of K ra . We assume that fl 
is bounded and that it has a C 2 -boundary. (The results hold for less restrictive 
assumptions.) We will use the index summation convention for repeated indices 
and subscripted comma followed by an index will indicate partial differentiation 
with respect to the corresponding variable. 

A vector field on Q is interpreted physically as a virtual velocity field on the 
body or alternatively as a field of virtual infinitesimal displacements. A rigid field 
in W is a vector field of the form 

w(x) = a + b x x, a,b G M 3 . 

Clearly, rigid fields may be restricted to subsets of R 3 . We denote the space of rigid 
fields by 2% and it is a 6-dimensional vector space. 

The following definitions and results concerning LZ3-fields are due to Temam 
and Strang |121 11UI ITT] . Given an integrable vector field w on f2, we consider the 
corresponding stretching (linear strain) field e(w) defined by 

e{w) lm = \{w^ m + w m .i), (2.1) 
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where comma implies distributional derivative relative to the corresponding spatial 
coordinate. The integrable vector field w is of integrable stretching, or w G LD(ft), 
if the components of the corresponding stretching are also integrable over ft. On 
the vector space LD(ft) of integrable stretchings it is natural to use the norm 

HI = \\M\ld = Nil + IkMlli, (2-2) 

where |j-|| p indicates the L p -norm. With this norm, LD(ft) is a Banach space and 
we have a continuous and linear 

s: LD(f2) — ► L x (.f2,R 6 ). (2.3) 

A basic theorem whose classical version is due to Liouville (see ^2 PP- 18-19]) 
states: 

Proposition 2.1. Kernel(e) = 

Let W be a Banach space of velocity fields. In the discussion below W will be 
either LD(ft) or L 1 (9i7, R 3 ). We refer to an element x € as a distortion. We 

have the natural projection mapping onto the quotient space 

tt:W — >W/M (2.4) 

and the induced norm in W/i? is given by 

||x||= inf||w||, or |[w]|| = inf||w + r||. (2.5) 

Proposition 2.2. For both W = L^&G.R 3 ) and W = LD(Q) there are continu- 
ous and linear projection mappings 

7r«:W— *&. (2.6) 

For w £ W, TTag(w) = a + b x x is given by 

' ' w, b = r L I / xx w \ , (2.7) 



\U\ 

u 

where 17 = dft for W = L^S^R 3 ), U = ft for W = LD(ft), \U\ is the Hausdorff 
measure of C/, and Ii m = JuixkXkfiim — XiX-m) is the moment of inertia of U. 

The space LD(ft) has the following properties (see 

Approximation: The restrictions of fields in C°°(.!?,R 3 ) to ft are dense in 
LD(ft). 

Extensions: There is a continuous linear extension operator E: LD(ft) — > 
LD(R 3 ). 

Regularity: If iu is any distribution on J? whose corresponding stretching is 
L\ then w G L 1 (r2,R 3 ). 

Trace mapping: There is a unique linear, surjective, continuous trace map- 
ping 

7 : LD(ft) — > L l (dft,R 3 ) (2.8) 

such that j(w\n) = w\dn for all continuous vector fields w defined on ft. 
Distortions: On LD(ft)/M, 

llxlU = IKx)lli (2-9) 
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is a norm that is equivalent to the quotient norm described above. Thus, 
there is a constant C(f2) (depending on fl only) such that for every w £ 
LD(Q) 

inf \\w + HI s: C(n)\\e(w)h. (2.10) 

The infimum is attainable, i.e., for each uu € LD(f2), there is a rigid motion 
Tq satisfying 

\\w + r \\ = inf ||u> + r|| < a(fl)||e(iu)||i (2.11) 

Equivalent norm: If p is a seminorm on LD(Q) such that p(r) — implies 
that r = for every rgf (so p is a norm on , then, 

p(w) + \\e(w)\\ 1 (2.12) 

is a norm on LD(f2) that is equivalent to the original norm. In particular, 
using the trace mapping one can take 

p(w) = ||7(w)||i,8fl, (2-13) 
so the following is a norm that is equivalent to original 12.2J1 

Nix =ll7Hlli,sn + l|eHHi. (2-14) 

Furthermore, one may use the projection ir^ : L 1 (i9f2, R 3 ) — > ^ as in Propo- 
sition [^21 and a norm ||-||^ to obtain the equivalent norm 

IMI© = lk# °t( w )IU + (2-15) 

Forces are regarded as elements of the dual spaces to the corresponding spaces 
of virtual velocities. So for a generic space of velocities W, a force F will be 
a member of W*. The evaluation F(w) is interpreted as virtual work, or virtual 
power, performed by the generalized force for the corresponding generalized velocity. 

In case the space W of velocities is an L p -space, 1 ^ p ^ oo, (e.g., L 1 (9J7,M 3 )) 
a force may be represented by an element of the corresponding L q space with 
q = p/{p — 1) through integration. We will use the same symbol for the force and 
its representing field. For example, for t G L 1 (9f2,]R 3 )* = L°°(9J?, R 3 ) we have 



t(w) = I t-w. (2.16) 

dQ 

A force F e W* acting on a body is equilibrated if F(r) — for all rff. An 
equilibrated force F is of the form F = ir*(F ) for some F £ (W/&)* . Furthermore, 
7r* is norm preserving, i.e., 

h*(Fo)\\ = HFoll, (2.17) 
so one can usually identify an equilibrated F with Fq. 

2.2. Generalized stress concentration factors and norms of trace map- 
pings. The central mathematical object that we find suitable for formulating the 
continuum mechanics problem, particularly, those notions related to stress concen- 
tration is (LD(Q)/3$)* — the dual to the space of LZ?-distortions. Specifically, as 
described below and in further detail in , elements of this space may be repre- 
sented by essentially bounded stress fields and on the other hand, the dual of the 
trace mapping associates an element of (LD(fl)/3$)* with any equilibrated bound- 
ary traction field. 
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Consider the composite mapping 71-07: LD(Q) — > L 1 (9J7, R 3 )/^. It is noted 
that for any r £ 7(7*) is a rigid motion on df2, hence, ir ° ^{w + r) = tt o "/(w) + 
7T o 7(r) = 7T o 7(1/;). Thus, we have a well defined mapping 

7/»: LD{0)im — ► L x {dfl, M?)/M, (2.18) 

given by j/&(x) — tt 7(10), for some w € x- Clearly, 

TT o 7 = (7/^) o TT. (2.19) 

In addition, 

||7(w) + r|| 1 = || 7 (w + r)||i || 7 || \\w + r\\ . (2.20) 

Hence, 

||[ 7 N]||= mf||7H+r||i<|| 7 || inf|k + r|| = ||7||||HII (2-21) 

r£3z r£3€ 

and we conclude that ^fjSt is indeed bounded and 

Win < IItII- (2-22) 

The dual mapping 7* : L°°(dQ 1 R 3 ) — > (LD(fi))* may now be applied to traction 
fields and (7/^)* : [L l (dfi, E 3 )/^)* -> (LD(Q)/M)* may be applied to equilibrated 
boundary traction fields to give LZ?-forces and equilibrated LZ?-forces respectively. 
Clearly, 

7T*o(7/^)* =7*0 7T*. (2.23) 

Proposition 2.3. The mappings 7* and (7/S?)* are injective. 

Proof. The mapping 7* is injective because 7 is continuous and surjective. As the 
quotient space projection 7r: L (e?J?,R 3 ) — * L (dQ,M. 3 )/M is also continuous and 
surjective, the same argument applies to j/M. □ 

Remark 2.4. Henceforth, we will use the equivalent norm ||-|| e as in l|2.9|) on 
LD(Q)/M. We will also use the equivalent norm 1 1 - 1 1 ^b> 011 LD(S7) as in (|2.15|) . This 
implies that the quotient norm on LD(ff)/& is actually equal (not only equivalent) 
to the norm induced by the strain. That is, 

HMIIHIMIIsHkHlli. (2.24) 

It is noted that the last equality, to be used below frequently, is independent of the 
choice of a particular projection ir^ : LD(f2) — > M. 

Proposition 2.5. Any T £ (LD(Sl)/&)* is represented by some symmetric stress 
field a e L°°(^,R 6 )*, in the form 

T=e*(a), (2.25) 
where e* : L°°(,G,R 6 ) -> (LD(SI)/^)* is the dual mapping to 

e:LD(n)/M — >L 1 (/2,R 6 ). (2.26) 
In addition, for the dual norm ||-|| e on (LD(f2)/&)* , we have 

inf JlalU, (2.27) 

it, T=e* (<r) 

and the infimum is attained for some a £ L°°(J7, R 6 ), i.e., 

II^THHU (2.28) 
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Proof. Using the duality L l (Q,R 6 )* = L°°(f2,R 6 ), the assertion follows from the 
fact that e: LD(Q")/M — > L 1 (/?,R 6 ) is a norm-preserving (by our choice of norm as 
in Remark 12. 4f> . linear injection and using the Hahn-Banach theorem (see [5] for 
the details). □ 

Corollary 2.6. Let t G L°°(df2, R 3 ) be any equilibrated traction field so there is 
a t G (L 1 (9J7, K 3 )/^ 1 )* such that t = 7T*(to). Then, there exists some stress field 
cr G L°°(f2,R 6 ) such that 

(7/«)*(to)=e*(a-), and 7* (t) = n* o e* (a) . (2.29) 

In addition, 

|| ( 7 /#)* (to) || = inf Halloo, (2.30) 

e»(<r)=(7/.^)*(*o) 

and 

Il7*(*)ll = M Iklloo- (2.31) 

7T* o £ *(cr)=7*(t) 

Proof. We make repetitive use of Ij2.23|l . I|2.17(l and Proposition H2.5JI . For example, 

H7*(*)ll = ll7*°7r*(to)|| 

= ||7r*o( 7 /«)*(t )|| 
= \\HW(t Q )\\ 

= mf HctHoo (2.32) 

it* ° £*(ct)=7t* o ( 7 /S?)*(t ) v ' 

= inf || cr|| oo 

it* o (cr)— 7* o 7T*(to) 

= mf Iklloo- 

7r* o e*(cr)=7*(t) 

□ 

Remark 2.7. The conditions i|2.29[) are equivalent to the principle of virtual work — 
a weak form of the equations of equilibrium — of continuum mechanics (as it is 
assumed throughout that the body forces vanish). For example, 7 *(t) = it* ° £*(cr), 
implies 

7 *(t)(w) = (7T*oe*)(a)(w), (2.33) 
so that for any w G LD{Q), 

t(~f(w)) = a(e o tt(w)) = <r(e(w)). (2.34) 
Hence, for a vector field «; that is the restriction of a differentiable field on J?, 



UWi = J aij£(w)ij. (2.35) 
do n 
Theorem 2.8. Let, 

||7/^|| = sup jj— J, , (2.36) 

x &LD{Q)/m 11X11 

be the norm of the trace mapping for distortions. Then, the stress concentration 
factor K for the boundary traction problem, satisfies 

K=hm\\- ( 2 -37) 
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Specifically, 



K = sup (2.38) 

»GC»(77) £ H l 



Proof. We have the standard 

|| 7 /^|| = ||(7/^)1- (2-39) 

However, 

IIWII=»pKMC (2.40, 

to ll*o|| 

sup j— inf Ikllool (2.41) 

toG((017,R 3 )/S?)* I 11*0 II £*(< :r ) = (7/^)*(to) J 

= sup {-L_ inf IHU}, (2.42) 

tGlmagcTT* I Plloo ?r* ° £* (<x)=7* (t) J 

where we used Proposition (|2.5(l and Corollary l|2.6(l . The condition 7r* o£*(cr) = 
7*(i) is equivalent to the condition that the stress field a is in equilibrium with t, 
i.e., <7 S by Remark 1)2. 7fl . The expression l|2.38[l simply uses the definition of 
\\j/S?\\ and the fact that C*°°(77) is dense in LD{Q). □ 

2.3. The scalar case and the trace mapping on the Sobolev space Wl(Q). 
The previous discussion is simplified considerably if we consider scalar fields tp € 
Wl(Q) instead of the vector fields w £ LD(S1). In this case the boundary data is 
also a scalar field, the analog of a stress is a vector field, and the vector space S? of 
rigid motions is replaced by the real numbers. 

Consider a sourceless vector field a on fi that satisfies boundary conditions for 
its boundary flux, i.e., 

a i;i = 0, in Q, (2.43) 
am = t, on df2, (2.44) 

for some given essentially bounded t : dfi — > !L Physically, a may be thought of as 
a material flow field (say for an incompressible flow) for a given flux density t on 
the boundary. Alternatively, a may be thought of as a heat flow field where there 
are no heat sources in Q so t is the given heat flux on the boundary; or er may be 
the electric displacement field and t is the charge density on the boundary. 
The weak formulation of the problem is 



= J tip. (2.45) 
n an 

(The test function <p may be thought of as a potential field in the electrostatic 
example or as the reciprocal of the temperature in the heat transfer example.) 
The analog of the stress concentration factor is then 

v ts ■ f IMLjJ . , ess sup xer2 |cr(a;)| 
K = sup K t = sup inf Tj-j: ! — = sup mf 2 — 

t£L°°(dO) t£L°°(dn) aeS t ||*||oo,9fi tGL°= (9J7) CTeI; t eSS SUp yeaj7 |<(?/) | 

(2.46) 

We will refer to K as the generalized field concentration factor. Thus for example, 
for the interpretation of a as a flow field, for a given flux density t, K t will be the 
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smallest ratio between the maximal magnitude of the velocity and the maximal 
value of the given boundary flux. 

Thus, ip may be regarded as an element of the Sobolev space Wl(fl) so 

IMIw? = IMIi + l|Vy>||i. (2.47) 

On the Sobolev space, the trace mapping 

7: Wl{Q) -» L l {dO) (2.48) 

is well defined as expected (see pQ). Our assumption that there are no sources in fl 
implies that J gn t = which is equivalent to considering Wl((l)/M. and 7/IBL Thus, 
the analog of Theorem l|2.8(l will be for the scalar case 

k = \\'y/r\\, y.w£({2) — yL^on). (2.49) 

3. Bounds on the W^-Trace Operator 

3.1. The bounds obtained using normal vector fields. In this section we 
consider bounds for the trace operator 

j:W { (f2) — > L}(dfl). (3.1) 

In particular, we are looking for bounds A and B satisfying 

J\v\^aJ\Vv\+bJ\^\ (3.2) 

dfl fl Q 

for every ip £ Wl(fl). 

Let n be any C^-vector field on fl and ip the restriction to Q of a VF^-function 
defined in an open neighborhood of fl. Then, 

rkip,i = J (ttiVO.-i _ / n i,i^ ( 3 - 3 ) 
on q 

implies using the Gauss-Green theorem that 

riiViip = J xiiip,i + J xii,iip, (3.4) 
dn n n 

where v is the outwards pointing unit normal to dfl. 

Definition 3.1. The vector field n on fl will be referred to as a normal field if the 
following conditions hold. 

(i) n(y) = v(y) for all y £ dfl. 

(ii) \xi\{x) ^ 1 for all x £ fl, where here and in the rest of this section we use 
the Euclidean norm for elements of W 1 so Inl = \/n7n7. 



The existence of normal vector fields is discussed in some detail below. We 
will use Vi{fl) for the collection of all normal vector fields. For a normal field, 
Equation (|3.4|l assumes the form 

ip = Jriiip^ + J n^itp. (3.5) 
dfi n a 
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Given a C 1 mapping (p on fl, the distributional derivatives \ip\^ of its absolute value 
= are clearly integrable and hence, \ip\ is W\. Rewriting Equation (13.51) 

for tjj = \ip\ we obtain 

J\<p\ = J m\<p\,i + J n iti \<p\. (3.6) 
an n n 
We now estimate each of the integrals on the right hand side. 



(3.7) 



»iiM,t ^ / MIM,il 

n Q 

Q 

where we used Definition (|3.1I m) and ||(/?| j| = |</5 j| ^ |Vy|. Also 

tli,i|¥?| «S max{|n M |} / (3.8) 

It follows that 

< / |V^| +max{|n M (x)|} f \tp\. (3.9) 
J xen J 

an n n 

This inequality is clearly exact and recalling the definition of normal vector fields 

we have 

Theorem 3.2. Let fl be a bounded open set in R™ having a C 2 -boundary and set 



B(fl) — inf < max{|tij i(x)\} 
Then, the following exact inequality holds 



K^U^M- (3 - 10) 



IMkan = J M < J |Vp| + B(fl) J \<p\ = \\V<p\\i + B(fl)|M|i. (3.11) 
an n n 

3.2. Estimation using harmonic normal fields. We now consider the existence 
of normal vector fields. By the assumption that dfl is C 2 , v is a C 1 - vector field on 
dfl and we can extend it to a vector field n with |n(a;)| ^ 1 for all x G fl (see for 
example Theorem 3.6.2 in J3j). Furthermore, we can require that the extension is 
harmonic in the following sense. For a vector field n, we use An for the vector field 
(An)j = (%■),«■ The field n is harmonic in fl if 

An = 0, for all x G fl. (3.12) 

Theorem 3.3. Let fl be a bounded open subset of R™ such that dfl is C 2 . Then, 
there exists a unique normal vector field no G 9T(J?) which is harmonic. In addition, 
for the supremum of the divergence, V • no = noi,i) we have 

l! V " n olloo.7? = H V ' n o||oo,or?- (3.13) 

Proof. For any fixed j, we have a classical Dirichlet problem 

Ar\j = in fl, n^ = Vj on dfl. (3-14) 

Given our smoothness assumption on dfl, there is a unique solution noj to each 
such boundary value problem and we obtain the harmonic vector field n . 
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For | no | 2 = nojtioj we have, 

A(n 0j noi) = ("ojtxojO.M 

= (2n j,,noi),i (3.15) 

= 2noj,Mnoj + 2no Ji inoj i i. 

In the last line, the first term vanishes because rioj is harmonic and hence A(|rio| 2 ) ^ 
0. We conclude that |no| 2 is subharmonic in Q. By the maximum principle for 
subharmonic functions 

max|n (V)| 2 = max|n (y)| 2 = max>(?/)| 2 = 1. (3.16) 

Thus, in addition to the boundary conditions, no satisfies the condition (|3.1I m). 
and so no € 0T(J?). 

Next, we note that for the harmonic no, 

A(V • no) = (noj,j),u 

= ( n oj,u),j (3-17) 
= 0. 

Thus, V • no is also harmonic in S7 so by the maximum principle Equation (|3.13|) 
holds. □ 

It turns out that the harmonic vector field of Theorem (|3.3() plays an important 
role in the computation of B(Q), the second constant of the sobolev Wl(fl) trace 
inequality. Continuing to use no for the unique harmonic normal vector field we 
have 

Theorem 3.4. The Sobolev constant B(S1) is given by 

B(f2)= inf ||V-n|| oo7 j = ||V-n ||cx 3 ,9fi. (3.18) 

Proof. Let (n m ), m € N, be a sequence of normal vector fields such that 

lim ||V • nJI = B(Q). (3.19) 

m — >oo 

Using normal tubular neighborhoods (e.g., jS] p. 110]), there is a 8 > such that 
we can parameterize an open neighborhood V of d£2 in ft by 

(y,z)€dOx[0,S), (3.20) 

where for each x € V , y(x) is a unique point on the boundary such that x is on 
the line through y which is normal to the boundary, and z is the distance to the 
boundary along the normal line where x is situated. For each m, let V m be the 
open set 

V m = {x e T2 : y(x) E 3D, z(x) < S/m} , (3.21) 

and let fi m = Q — V m so 

dQ m = {x e T2 : y(x) € <9i7, z{x) = 5/m} . (3.22) 

Now for each m we construct the harmonic lifting (cf. ^ p. 24]) n m of n m as 
follows. Let nom be the solution of the Dirichlet problem on Q m with the boundary 
conditions no m (a;) = n m (x) for x € df2 m - Set 

, , tTn m (a^), for x G i7 m , 
ri m (x) = { 0m //' _ m ' (3.23) 

n m (a;), for x € V m . 
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By the maximum principle, 

IKrflloo.U < ll n »»lloo,77 and ll V ' "mlloo.U < H V ' n ™Hoo,7?> ( 3 - 24 ) 

so 

Jm j ||V-n Tn || 00i7 j = J B(fl). (3.25) 

In addition, as the various lltwll^ are bounded by 1, the same applies to n m i 
and the sequence n m i is uniformly bounded. Thus, (using a standard normal family 
argument) by Ascoli's theorem, it has a subsequence that converges uniformly to a 
limit continuous normal field n. On any compact subset of fl this gives a uniformly 
convergent sequence of harmonic functions whose limit is then also harmonic. Thus, 
n is harmonic. Also, the limit n satisfies the conditions of Definition (|3.1(1 and so it 
is a normal vector field. Finally, by the uniqueness of the solution to the Dirichlct 
problem n = no . □ 

Remark 3.5. Within the framework of the Ai?-program in geometric analysis, 
M. Motron proves in using different methods, the following two theorems for 
bounds on the trace mapping on Wl(fi). 

(1) For any e > 0, there exists a B E such that for any tp £ Wl(f2), 



+ /IV^+Se /M- (3.26) 

an a n 

(2) Assuming that Q is a connected bounded open subset of R™ whose boundary 
is piecewise C\ there exists A > such that for any ip £ Wl{Q), 

\<p\^A /|V^| + M / M . (3.27) 



an n a 

We note that even if we had the solution to the v4B-program, we would not have 
a concrete bound on ||7||. What we sought were simultaneous bounds on both A 
and B. 

In addition, for a normal vector field 

sup|V • n||J?| ^ /V-n= fn-v = \df2\, (3.28) 
xen J J 



i).Q 



so 



\en\ 

ni5t" " " llo °'" Ti? 



^lllV-nlU^^^. (3.29) 



However, the inequality is not exact and equality is not attainable even for the 
harmonic normal vector field. (Think of two circles connected by a narrow neck of 
width t. Across the narrow neck, V • n has to be of order 1/t in order to satisfy 
the boundary conditions. The values of |4?| and \df2\ are not significantly different 
from those of the two circles.) Thus, Motron's bound B is smaller than the value 
obtained here. On the other hand, the bounds A and B we obtain are exact in the 
sense that you cannot lower B without increasing A. 
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4. Norms of Symmetric Tensors 

As noted earlier, the value of the stress concentration factor depends on the norm 
we choose to use on the space of stress matrices. Thus, for the sake of completeness, 
we review below some elementary properties of norms of symmetric matrices on R". 
We will denote the norm of a matrix T by \T\ (reserving ||-|| for norms on function 
spaces). 

4.1. Operator norms. In general, for a linear mapping T: V — > U between 
normed spaces, the operator norm of T is defined by 

|T| = sup^l, v?0. (4.1) 

v \V\ 

The operator p-norm, |T| op , 1 ^ p ^ oo, on the space of matrices is defined as the 
operator norm for the case where the p-norm is used on both V = R m and U = 1™. 
By the compactness of the unit ball in R m , the supremum is attainable and 

|T| op = max|7»| p . (4.2) 

\v\ — 1 

In case T: W — > W is a linear transformation defined on the inner product 
space R™ equipped with the Euclidean 2-norm, |T| Q 2 may be calculated by 

\T\ 2= sup v,v'^0. (4.3) 

v,v'EW \V\\V'\ 

The following relations hold for symmetric matrices on R™. 

|T| ol = |T| 000 =max^|Ti j |, (4.4) 

j 

|T| o2 =max{|A 1 |,...,|A„|}, (4.5) 

where Ai, . . . , A„ are the (real) eigenvalues of T. The norm \-\ o2 is usually referred 
to as the spectral radius norm. 

4.2. Vector norms. We will also regard symmetric matrices as vectors in R™(™+ 1 )/ 2 
and use the p-norm for them. Thus, 



i,3 



In particular, 



\T\i = J2\ T v\> ( 4 - 7 ) 

|T| 00 = sup|T ij |, (4.8) 

i,3 

|T| 2 = v/7p>- (E A ?) 1/2 - (4 - 9) 

We recall that the Frobcnius norm of a matrix is |T|i? = {TijTij) 1 / 2 and it is 
identical to |-| 2 for symmetric matrices. 
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4.3. Dual norms. Being a finite dimensional space we may identify the space 
of symmetric matrices R n ( n + 1 )/ 2 with its dual space. Thus, we may regard any 
symmetric matrix T as a linear functional so that T{S) = TijSji and assign to it 
the dual norm |T| P » 

|T| p ,=sup^pl, (4.10) 
s Pip 

where we have the usual |-| p * = |-| 9 , for q=p/(p— 1). 

Dual norms may be used also for the operator norms. In particular, we note 
that 

|T| o2 .=5>|. (4.11) 

i 

In closing this short review, it is noted that the norms containing the index 2 
are associated with the Euclidean norm for vectors and may be expressed in terms 
of the eigenvalues. These norms are invariant under orthogonal transformations of 
coordinates. 

4.4. The equivalence constants. Since all the norms listed above are equivalent, 
for each pair of norms |*| a cmd |*|^, there is & finite positive number 

tf 6 » = 8up{^. (4.12) 
T \T\ b 

In particular, the following exact relations hold: 

V n \ T \ol \T\o2 \T\oc 

1 < < n, 1 < < ^. (4.14) 

I 1 \o2 K \2 

Thus, for example, K % 2 = 1. 

In the sequel, we will use |cr(x)| to denote the norm of the value of a symmetric 
tensor a at x € O using a generic norm on the space of matrices. 

5. Optimal Boundary Conditions for Stresses 

Consider the following problem: given a unit vector v in a Euclidean 3-dimensional 
space and a unit vector t, find a symmetric matrix a such that 
(i) <j(v) = t — the compatibility condition; 

(ii) \<t\ = inf{|T|, T{v) = t,T = T T }, i.e., a is the optimal symmetric matrix 
that satisfies condition (i). 

The problem has an obvious mechanical interpretation. If v denotes the normal 
to the boundary at some given point, and t denotes the value of the surface traction 
field at that point, then, a(v) = t is the boundary condition for the stress field a. 
Thus, a matrix a satisfying the conditions above is the optimal stress matrix that 
will satisfy the boundary condition. Obviously, the normalization condition on t 
causes no loss of generality. 

Let t n = (t ■ v)v — v® v(t) be the normal component of t and let t t — v x (t x v) 
be the tangent component of t. Thus, denoting the angle between v and t by 6, 
\t n \ = cos# and \t t \ = sin#. We choose a basis {fj} where /i = v, fi is a unit 
vector in the direction of tt and /jj completes the other two to form a right-hand 
oriented orthonormal basis. In this basis, the matrix of a that satisfies the condition 
a(u) = t has to satisfy <7n = cos#, cri 2 = <j 2 i = sin^, and a\3 = 031 = 0. The 
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rest of the components cannot be determined by the compatibility condition above 
and should be determined by the requirement for minimal norm of a. (In the case 
where v and t are parallel, one can take any orthonormal basis containing v.) 

5.1. Optimal boundary conditions relative to the |- loo-norm. We wish to 
regard a as an element of the dual space of symmetric matrices. Then, using 
the basis fi as above, the compatibility condition implies that we have a linear 
functional do defined on the subspace V of symmetric matrices containing elements 
of the form 

en £12 

[e]= ei2 . (5.1) 


The functional ctq acts on elements of V by 



a o([e]) = £n cos9 + 2£i2 sin 9. 



Thus, on this subspace 

sup{|(7 (£) 

The extension 



;e G V, |e|x = 1} = max{|cos0|, |sin#|}. 



(5.2) 



(5.3) 



(5.4) 



cos 8 sin 8 
sin 8 


of do to the space of all symmetric matrices, has the same norm and as such it 
provides the optimal boundary condition. It is noted that while the development 
depends on the basis chosen, the optimal norm depends only on the angle 8, an 
invariant quantity. 

5.2. Optimal boundary conditions relative to the | |2-norm. If one uses 
M = 1 17 1 2 induced by the inner product in the space of symmetric matrices as in 
the previous section, the optimal matrix can be obtained by orthogonality condi- 
tions of the optimal stress to basis vectors for the matrices that correspond to the 
undetermined components of the stress. This implies that all the undetermined 
components should vanish. Thus, the optimal stress is given in the {/i}-basis by 
Equation (|5.4J) also, and 

\a\ 2 = Vl + sin 2 6. (5.5) 

While the optimal a for the | -loo-norm depended on the basis {fi} chosen, the 
construction here is rotation- invariant. 

5.3. The case where t — e^. We now consider the special case where t = e k , where 
e k , k — 1, 2, 3 is a base vector. In this case, cos 8 — v ■ e k = vu and sin 8 = a/1 — v\. 
Using 

cr = cos 8v ® v + sin 8{v ® f 2 + f 2 <8 v), (5.6) 
we may write the matrix for a relative to the {ei}-basis. We first note that 

e k -cosdv e k -v k v 



h 



sm ( 



0" 



(5.7) 



thus, 



a = v k v ® v + v ® (e*; — v k v) + (e k — v k v) % v. (5.8) 
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Rearranging the terms we conclude that for the case t — eu, the optimal boundary 
conditions for the stress are 

a = — Vkv (g) v + {v ® ejt + efc ® v), (5.9) 

and the optimal values are 



M 2 = V 2 -^ 

| a loo — max -||^A;|, \J 1 — v 2 j- , for the natural basis {/j}. 

5.4. Example: the 2-dimensional case for the |-| 2-spectral radius norm. 

Using the same notation as above and using the coordinate system in the plane 
where the x and y axes are along f\ and /2, respectively, we are looking for a 2 x 2 
symmetric matrix that will satisfy the condition t — a(y) of least spectral radius, 
i.e., minimizes max{|Ai|, | As j } - The condition a(v) = t implies that a xx — cos# 
and a xy = sin 9. Thus, to determine a completely, one has to determine the single 
number a yy . 

In two dimensions we have the explicit expression for the eigenvalues as 




*1,2= VV ±\,[ " "" ) +r,,~ (.-,101 

Setting 



a= cos9 = ^^ z= ^ ; so al y =sin 2 e = A-a 2 , (5.11) 



we have 



Ai, 2 = a + z ± y/l - 3a 2 - 2az + z 2 = a + z ± y/{z - a) 2 + 1 - 4a 2 . (5.12) 
Minimizing | Ai,2 1 is like minimizing \\ 2 so we differentiate with respect to z. Thus, 
^1.2 d\i j2 

— — ZAi 9 

= 2Ai.2 



V(z - a) 2 + 1 - 4a 2 y 
that vanishes identically only if 4a 2 = 1. Hence, for extremum, 6 = 0. In general, 

1 -4a 2 = 1 -cos 2 6> > 0, (5.13) 

which implies that 

-ff (5-14) 

az 

has the same sign as the eigenvalue. It follows that z / can only make | A1.2 1 
larger and the infimum is attained for z = a yy = 0. 

Remark 5.1. We note that the problem of optimal boundary condition for the stress 
is a generalization of the requirement in Definition (|3.1f i)) for the boundary value 
of a normal vector field with the difference that now we consider matrices rather 
then vectors. Indeed, n = v gives the smallest value for | n| 2 , the Euclidean norm of 
the vector field n, such that n • v = n(v) = 1. 
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5.5. Worst case optimal boundary conditions. By worst case optimal bound- 
ary conditions we refer to the inclination of t relative to the normal v for which the 
norm of the optimal stress attains a maximal value. That is, we are looking for 

D|| = sup jinf |cr| : a{v) = t, \t\ = 1 j . (5.15) 

The number Du (the subscript || indicating the particular norm chosen) and the 
corresponding a and t depend only on the choice of norms. For example, for the 
|-|2-norm, D2 = V2 is attained for any unit traction t perpendicular to v. For the 
|-|oo-norm, relative to the boundary natural basis, Dqo = 1 is attained for cither 
traction that is parallel to v or traction that is perpendicular to it. 

5.6. The spaces 9T fe (J2) and e^-optimal tensor fields. For k = 1,2,3 we now 
choose a symmetric tensor field T k — (T* ) on dfi that satisfies the following con- 
ditions. 

(i) If e k denotes the fe-th base vector in R 3 , then, T k {v) = e k so T k -v^ — S k . 

(ii) Clearly, at each point y on the boundary there is a collection of matrices 
that satisfy condition (i) above. We choose T k (y) to be a symmetric matrix 
that satisfies condition (i) and that has the least norm (of our choice on 
the space of matrices). Thus, 

\T k (y)\=ig{\S k \:S k (u(y)) = e k }, (5.16) 

where the infimum is taken over all symmetric matrices. Thus, in the 
terminology of the preceding subsections, T k is the optimal stress for the 
efc as boundary traction and the discussion of Subsection i|5.3[l applies. 

Clearly, the fields T k depend on Q only. They depend continuously on v so by 
the assumption that dfl is C 2 they are continuous. Consider now the worst value 
of \T k (y)\ on the boundary, i.e., 

sup {\T k (y)\}= sup \w£{\S k \:S k (v(y)) = e k }\. (5.17) 

yean y ean I s k ) 

As dfi is assumed to be smooth, any angle between v and any fixed vector is 
attained on the boundary, hence, the worst case optimal boundary conditions are 
attained on the boundary always. Thus, 

sup {|T fc (y)|}= sup {mf{\S k \:S k (v(y))=e k }\=D ll (5.18) 

and depends only on the choice of norm. 

We use the notation T k for a tensor field on the boundary satisfying the two 
conditions above. Using Whitney's extension theorem (cf. ^3 Theorem 3.6.2]), T k 
can be extended to symmetric differentiable tensor fields cr k on Q that satisfy the 
following condition 

sup {|«t*(x)|} = sup {\T k {y)\} = D||. (5.19) 

We denote the class of symmetric C^-tensor fields u k on Q that satisfy the 
boundary condition cr (y) = T k (y), y £ df2, and condition l|5.19|l above by 9T fc (i7). 
We will refer to a tensor field cr k S 9t fe (J?) as an ek-optimal tensor field. 

It is quite clear that the foregoing discussion applies in the continuum mechanics 
context to the optimal stress field for the normalized boundary traction t = e&. 
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If a k 6 7t k (f2), we have for the stress concentration factor and optimal stress 
concentration factor, 

K ehttF >=K eh =D p (5.20) 
Although to total force on the body is not equilibrated, Dm, may serve as a bound. 

6. Bounds of the LD(J?)-Trace Operator 
Theorem 6.1. Let the constants A(f2) and B(fl) be given by 

A(/2) = 3D||, (6.1) 

B{Q) = Y)W ■ 4\\°o,aa, (6.2) 

k 

where (Tg is the solution of the Dirichlet problem Acr^ = 0, in J?, er*j = T* , on df2. 
Then, 

\\w\\i, 9 q <4(fl)||e(ti0||i + D(^)H|i, (6.3) 
for all u> € LD(f2), is an exact estimate. 

The following subsections present the proof. (We will use A and B for j4(J?) and 
B(Q), respectively, in order to simplify the notation.) 

6.1. The principle of virtual work. Let Q be an open region in R 3 having a 
smooth boundary, <j = (cy) a symmetric smooth tensor field on fi and w — (wi) 
an LD-vector field on fl. Then, 

Also, by the symmetry of a, 

PijWij = (TijEij. (6.5) 

Thus, we may write 

" j^,^ ■ ja,,,r,. (6.6) 
n n n 

and using the Green-Gauss theorem on the first term on the right-hand side we 
obtain 

1 



a ij £ ij = J "iig^.i + W 3,i) = J 0-ijWiPj+ J (TijjWi, (6.7) 

q n an n 

where v is the unit outward pointing normal. We will refer to the identity above 
as the principle of virtual work. 

6.2. The bounds. We now write the principle of virtual work (j6.7|l for the vector 
field (w) — {\wi\). (We reserve the notation \v\ for the norm of the vector v in a 
finite dimensional space.) Thus, 

Oij\{\ w i\,i + \ w j\,i) = J Oy-Kkj + J <rij,j\u>i\- (6-8) 
n an n 

Let a k satisfy a k (v) = ek on dQ so cr^Vj = S k . Then, the identity above assumes 
the form 

hk|= f\{\<j + \wj\,i)v%-j°h\™i\- (6-9) 
an n n 
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Consider the integrand 



1 



integrand = ^(\wi\j + \wj\,i)a^ (6.10) 

of the first integral on the right. As the expression is invariant under orthogonal 
transformations, it may be evaluated in the principle coordinate system of the 
matrix (|u>i|,j + \wj\ t i)/2 where the off-diagonal elements vanish. Thus, without 
loss of generality, we may write (we do not use the summation convention here) 

integrand = ^KLiof; 

maxjJofjiyVil.i 

i 

< W k \oo ^2,\ w iA (using \wi\j = sign(w i )w lji ) 

i 

<k fc LK«0li. 

where the equality is clearly attainable. Thus, 

/KK Jw k L\e(w)\ 1 + J \4,j\\ Wi \ 
an n n 

sup{|a fe (a;)L} f\e{w)\ 1 + sup_{\^ jtj (x)\} f ^Kl- 

As 

™v_\<J* jJ {x)\ = \\V-o k \\ 00 j S , (6.11) 

i,x Ef2 

we have 

\\Mi,an = J Kl < II^IL^IkHlli + IIV ■ ^H^HIi, (6-12) 

on 

where we use \\T\loo = HITIooHloo, and ||T||i = |||T|i||j,i, for the respective norms 
of a tensor field T. 

Adding this equation for k = 1, 2, 3 we obtain for the Z^-norm of the restriction 
of w to the boundary the following bound 

\\w\\ 1M2 < EII^IUt^HIIi + E !l v • ^IL,7?lklli- (6.i3) 

k k 

Clearly, for 

A = mf|x;ik fc lloc,fl : = (6.14) 

and 

B = 5? |EII V • ^Hoo.TJ : = e fe | . (6.15) 
this bound is the tightest and we have 

IIHkafl < A \H W )\W + B \\M\u for a11 weLD(H). (6.16) 
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Since in general 

sup|cr fc (x)| QO > sup \o k {y)\oo 
xel2 yedo 

^ sup (T^^loo, T k optimal as in fiTty 

yedn 

= D \\ 

and equality holds for a k € we conclude that A is attained for fields a k G 

^(ft) and A = 3D||. This proves the first part (Equation iJJQJ) of Theorem IO 

6.3. Estimating B. The procedure we use is completely analogous to that of Sub- 
section 0Ol an d the proofs of theorem !3.3l and Theorem 13. 41 For a tensor field a we 
use Act for the Laplacian Act^ = Uij.u- We say that a is harmonic if Ac = 0. 

Proposition 6.2. There is a unique harmonic tensor field a k 6 yi k (£2), i.e., a k is 
efc-optimal. For the harmonic e^-optimal a k , we have 

\\V-<T k \\ 00 ,a = \\V-o- k \\ 00 , d n. (6.17) 

Proof. Consider the Dirichlet problem 

Aa k = 0, in Q, a k = T fc , on df2. (6.18) 

The existence and uniqueness are standard. Let <r k be harmonic, then, for each 
component 4 (no sum on repeated indices) 



a((4) 2 ) = (44), 



2(4 /4i) ; 

' ' (6.19) 



> 0. 

Thus, 4 i s subharmonic in £1. By the maximum principle for subharmonic func- 
tions 



max (4 (x)) 2 = max {a% (y)) 2 (6.20) 
and so the analogous property holds for \o~ k A. Thus, using the boundary conditions 



lk & lloo,77 = ma*J4(*)l = ■ max 1^.(^)1 = . max jT^y)] = D x . (6.21) 

i,j,x£f2 l)3>y£ot& i.j,yEuiJ 

Hence, the solution is also in 9t fc (J2). Finally, 

A(V-a fc ) = (4 )i ),H 

= (4,u)j ( 6 - 22 ) 

= 0, 

so by the maximum principle for the components of V • a k , (|6.17() holds. □ 

Proof of the second part of Theorem \6 . 1\ (Equation \6. ty) ). Clearly, as the three o~ k 
fields are independent, we should look for 



B k 



^ f {ll V ' ^HooJJ :ff*M=e fc }. (6.23) 
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Thus, let (cr^), m € N, be a sequence of e^-optimal tensor fields such that 

lim ||V • a* || = B k . (6.24) 

m — >oo 

The sets V m and f2 m may be constructed just as in the proof of Theorem 13.41 Let 
<7o m be the solution of the Dirichlet problem in fi m with the boundary condition 
a 0m( x ) = a rn( x )i f° r x € df2 m , and define the harmonic lifting accordingly as 



0omO)> for x e n, 

<Tm(x), for x g V 



w) = rr:~ zzZtt (6-25) 



By the maximum principle 



<ij\L,-n ^ \\<ij\L& and l|V-<ILj7< IIV-^IL,^, (6.26) 
lim j|V •<|| 0O o = S fc . (6.27) 



m — >oo 



We apply the normal family argument and uniqueness of solution as in the proof 
of Theorem l3.4l to obtain 

Bfc = ||V-(rg||oo,sn, (6-28) 
where Cq is the solution of the Dirichlet problem (|6.18() . Equation (|6.2() now follows 
from (|6"T5|) . □ 

7. Discussion 

We may describe the foregoing analysis in the mechanical context. For a given 
boundary traction field, we constructed the optimal boundary condition for the 
stress field. In Section 03 the case where the traction vector was a unit vector was 
considered, but as the relation between stress and traction is linear, this causes no 
loss of generality. Thus, one can assign the optimal boundary condition for the 
stress field for any given boundary traction field. 

Next, one can solve the Laplace equation for each of the stress components. 
Unlike the usual case of continuum mechanics, we have a unique solution without 
imposing constitutive relations and the equilibrium equations are not satisfied. For 
the harmonic solution of the boundary value problem, the maximal stresses occur on 
the boundary and these stresses are the smallest that satisfy the traction boundary 
conditions. Eg uat ion 16.21 and its proof indicate that the maximal value of V • a is 
the smallest possible. In light of the usual equilibrium equations V • a + b = of 
continuum mechanics (b being the body force field), we can interpret the field — V-cr 
as additional body forces one has to supply for the equilibrium condition to hold. 
Thus, for equilibrium, the harmonic stress field is associated with an additional 
body force field whose maximum is the least (and is attained on the boundary). It 
is noted that the total of the traction field, J gn t, was not required to vanish so it 
is not possible for equilibrium to hold. With the foregoing limitations in mind, the 
harmonic stress field solves the problem of optimal stress field for a given traction. 

Next, we note that the generalized stress concentration factor may be described 
as the largest optimal stress concentration factor when we can vary the boundary 
traction fields while keeping their maximal value on the boundary to ||t||oo,9fi = 1- 
Thus, in the analysis all the components of the traction fields are set to be 1 
everywhere. Again, this precludes equilibrium as the total force on the body in 
each direction is equal to the area of the boundary. The mathematical analog of 
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this limitation is that we obtain bounds on 7 and not j/M. However, the bound on 
H7II gives a bound 1 1 7/^ 1 1 because < ||7|| (Equation (J223). 
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